Measurement of the two neutrino double beta decay half-life of Zr-96 and search 
for associated neutrinoless processes with the NEMO-3 detector 



J. Argyriades^', R. Arnold'', C. Augier^, J. Baker*^, A.S. Barabash*^, A. Basharina-Freshville', M. Bongrand'', 
G. Broudin-BayS'^, V. Brudanin'', A.J. Caffrey'^, A. Chapon", E. Chauveau^'^, Z. Daraktchieva', D. Durand", 
V. Egorov'*, N. Fatemi-Ghomij, R. Flack', B. Guillon", Ph. Huberts-^, S. Jullian^, M. Kauer*'*, S. King', 
A. Klimenko'', O. Kochetov'', S.I. Konovalov'^, V. Kovalenko'', D. Lalanne^, T. Lamhamdi'^, K. Lang', Y. Lemiere", 
C. Longuemare", G. Lutter^-^, F. Mamedov™, Ch. Marquet^ f, J. Martin-Albo°, F. Mauger", A. Nachab^-f, I. NastevaJ, 

I. Nemchenok'', CH. NguyenS-f ", F. NovaP, R Novella", H. Ohsumii, R.B. Pahlka', R Perrots-^, F. PiquemaF-f, 
J.L. Reyss', J.S. RicoF '^, R. Saakyan', X. Sarazin% Yu. Shitov'', L. Simard% F Simkovic\ A. Smolnikov'', S. SnowJ, 
S. Soldner-RemboldJ, I. Stekl"", J. Suhonen', C.S. Sutton", G. Szklarz% J. Thomas', V. Timkin\ V.I. Tretyak''\ 
V. Umatov^ L. Vala"", I. Vanyushin^ V. VasiUev', V. Vorober, Ts. Vyiov'' 

On 

(The NEMOS Collaboration) 
O 

^ "lAL. Universite Paris-Sud 11, CNRS/1N2P3, F-91405 Orsay, France 

Q ''IPHC, Universite de Strasbourg, CNRS/1N2P3, F-67037 Strasbourg, France 

D '^INL, Idaho National Laboratory, 83415 Idaho Falls, USA 

''INR, Institute of Nuclear Research, MSP 03680 Kyiv, Ukraine 
'ITEP, Institute of Theoretical and Experimental Physics, 1 1 7259 Moscow, Russia 
fCNRS/IN2P3, Centre d' Etudes Nucleaires de Bordeaux Gradignan, UMR5797, F-33175 Gradignan, France 
^Universite Bordeaux, CENBG, UMR 5797, F-33175 Gradignan, France 

1 I ''JINR, Joint Institute for Nuclear Research, I4I980 Dubna, Russia 
^ 'University College London, WCIE 6BT London, United Kingdom 

'University of Manchester, Ml 3 9PL Manchester, United Kingdom 
' I ^USMBA, Universite Sidi Mohamed Ben Ahdelkih, 30000 Fes, Morocco 

^ ' University of Texas at Austin, 78712-0264 Austin, Texas, USA 

"'lEAP, Czech Technical University in Prague, CZ-12800 Prague, Czech Republic 
"LPC, ENSICAEN, Universite de Caen, CNRS/1N2P3, F-14032 Caen, France 
I I "IFIC, CSIC - Universitat de Valencia, Valencia, Spain 

^ Universitat Autdnoma Barcelona, Barcelona, Spain 
CO ''Saga University, Saga 840-8502, Japan 

'LSCE, CNRS F-91190 Gifsur- Yvette, France 
'FMFI, Commenius University, SK-842 48 Bratislava, Slovakia 
'Jyvdskyld University, 40351 Jyvdskyld, Finland 
\ ^ "MHC, Mount Holyoke College, 01075 South Hadley, Massachusetts, USA 

^Charles University in Prague, Faculty of Mathematics arul Physics, CZ-12116 Prague, Czech Republic 
• ''Hanoi University of Science, Hanoi, Vietnam 

o 

On 

O 

^ Abstract 

^ Using 9.4 g of '*Zr and 1221 days of data from the NEMO-3 detector corresponding to 0.03 1 kg-y, the obtained 2vpp 

^ decay half-life measurement is T^J'^ - [2.35 + O.I4(.?fflf) ± 0.16(.s->\s-f)] X lO'^yr. Different characteristics of the final 

state electrons have been studied, such as the energy sum, individual electron energy, and angular distribution. The 2v 
nuclear matrix element is extracted using the measured 2v/3p half-hfe and is M^" = 0.049 ± 0.002. A 90% CL Umit 
is set on the Ov/3/3 decay half-life of T'^J^ > 9.2 x 10^' yr corresponding to a limit on the effective Majorana neutrino 
mass of {tri/jis} < 7.2 - 19.5 eV. Limits on other mechanisms of OvySyS decay have also been set. 
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1. Introduction 



The recent observation of neutrino flavor oscillations and the resulting measurements of the neutrino mass squared 
differences ID have motivated renewed experimental efforts to measure the absolute neutrino mass. The fundamental 
Dirac or Majorana f2] nature of the neutrino also remains indeterminate. Neutrinoless double beta decay {Ov/3/3) is the 
only practical means of determining the nature of the neutrino and one of the most sensitive probes of its absolute mass. 
The mechanism most commonly discussed is shown in Figure[T]in which a light Majorana neutrino is exchanged |[3l . 
The half-life in this case is given by 

[rf;^]-' =G»''|M"''|2(m^^)2, (1) 

where G*''' is the precisely calculable phase-space factor (proportional to Q^p), M"'' is the nuclear matrix element and 
(m^^) is the effective Majorana neutrino mass. Other possible mechanisms for OvySyS include, for example, right-handed 
currents, Majoron emission and R-parity violating supersymmetry. In all mechanisms, the Qv/3/3 process violates lepton 
number and is a direct probe for physics beyond the Standard Model. 




Figure 1: Feynman diagram for Ov/}fl decay via light Majorana neutrino exchange. 



Measurement of the 2v/3/3 process is important because it is an irreducible background component to Ov mech- 
anisms. Double beta decay (2vy6y6) allows experimental determination of the NME for this process (M^''), which 
leads to the development of theoretical schemes for NME calculations for 2v/3/3 and Ov/3/3 |4, 5. 6|. The precision 
with which the lepton number violating parameter, such as {mp/j), can be measured depends crucially on knowledge 
of M^^ . Presented here are the results of observations of ^''Zr obtained with the NEMO-3 tracking plus calorimeter 
detector. 



2. NEMO-3 Experimental Apparatus 

A detailed description of the NEMO-3 detector and its performance can be found in Q, while the most salient 
properties are mentioned here. The detector is located in the Modane Underground Laboratory (LSM) 4800 me- 
ters water equivalent below ground and has been acquiring data since February 2003. It is a cylindrical detector 
(0 5 X 2.5 m) holding 10 kg of enriched isotopes. The tracking volume contains ~ 6000 drift cells operating in Geiger 
mode (Geiger cells) enclosed by ~ 2000 polystyrene scintillator blocks making up the calorimeter The detector is 
enclosed in a solenoid which generates a 25 Gauss magnetic field parallel to the Geiger cells. The transverse and 
longitudinal resolution of the tracker is 0.6 mm and 0.3 cm (cr) respectively. The calorimeter energy resolution and 
timing resolution is 14-17% (FWHM at 1 MeV) and 250 ps (cr at 1 MeV) respectively. 

The majority of the /3/3 isotope mass is "'"Mo but other isotopes include '^^Se, "^Cd, '^°Te, '^"Nd, ^''Zr, and "^^Ca. 
The experimental signature of Qvf3/3 is two electrons with the energy sum equaling the Qpp of the decay. '*Zr is of 
particular interest due to its high Qpp - 3350.0+3.5 keV which is greater than the decay energies of most contributing 
background sources, and the large phase-space factor which is proportional to The total mass of the enriched 
ZrOz is 22.0 g of which 9.4 + 0.2 g is ''''Zr Q. NEMO-3 results thus far are pubUshed in ll8ll9l[T0l[ni. 
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3. Event Topology and Particle Identification 

The NEMO-3 detector is capable of sophisticated particle identification and event topology reconstruction. Elec- 
trons and positrons produce signals in both the calorimeter and Geiger cell tracker, while photons only signal the 
calorimeter Due to the 25 Gauss magnetic field permeating the detector volume, the electron and positron discrimi- 
nation efliciency is 97% at 1 MeV. Alpha particles {a) are identified by the short distance (~ 20 cm) they travel before 
quenching in the gas volume of the Geiger cells. Crossing electrons (an electron crossing the whole tracker volume 
and source foil to mimic a /3/3 event) are identified by the time-of-flight information from the two signaled calorimeter 
blocks. The event topologies studied in this analysis include the single electron channel (le), the electron plus gamma 
channel (ey), and the two electron channel (ee). 

4. Backgrounds in NEMO-3 

Studies have been carried out to identify the activities of the contributing backgrounds to NEMO-3 fT2l. The 
backgrounds are categorized as "internal" or "external." Internal backgrounds include isotopes decaying from within 
the source foil miinicking a ySjS decay via M0ller scattering, /3 decay with internal conversion, or/? decay with Comp- 
ton scattering of the de-excitation photon. Each of the seven /3/3 isotopes being measured at NEMO-3 has specific 
dominant internal backgrounds. External backgrounds include all decays originating from outside the source foil but 
still mimic a /3j3 event signature via double Compton scattering, Compton plus M0ller scattering, or pair production. 
Charge identification via track curvature in the magnetic field tags pair production events. The two most detrimental 
contributers are ^'"'Bi and ^"^Tl with respective g/J values of 3.270 MeV and 4.992 MeV. NEMO-3 component and 
source foil activities were measured with a high purity germanium detector (HPGe) and were subject to a selection 
process to optimize radio-purity. 

4.1. Radon (^^^Rn) 

The first data acquisition period (Feb 2003 - Oct 2004) is referred to as Phase-I and had a relatively high 
(1200 mBq) level of radon in the tracking volume. Radon (^^^Rn) is particularly disruptive because it is a noble 
gas and its half-life of 3.82 days provides enough time to be outgassed from the surrounding rock and permeate the 
detector volume. Supporting evidence suggests iTfTI that a large fraction (87%) of a decay daughters are positively 
charged and are attracted to electrically negative and grounded surfaces. NEMO-3 data are consistent with the radon 
daughters being deposited on the surfaces of reflecting wrapping around the scintillators, the drift cell cathode wires 
and the source foils lfT2l . 

The second data acquisition period (Nov 2004 - Dec 2007) is referred to as Phase-II and began with the installation 
of a radon purification facility to inject a flow of pure air around the detector. The purification facility suppresses the 
radon concentration in the immediate proximity of the detector by a factor of ~ 1000. However, the outgassing of 
detector components releasing radon due to their internal contamination with the ^^**U- ^^^Ra chain leads to a smaller 
reduction factor inside the detector. The radon activity in the tracker volume decreased from 1 .2 Bq in Phase-I to 
0.2Bq in Phase-II. 

5. Data Analysis 

All background and signal events are simulated with DECAYO (T4\ which accurately reproduces energy and 
angular distributions of particles emitted in radioactive decays including 2vy6y6 and theoretical Qv/3/3 mechanisms. All 
generated particles are propagated through a full GEANT-3.21 ifTSl description of the detector. The simulated Monte 
Carlo (MC) events are in the same format as the raw data from the NEMO-3 detector and both MC and real data are 
reconstructed with the same software package. 
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5.1. Background Identification 

One can measure the activities of the various background isotopes by the event topologies and kinematics deter- 
mined by the selection criteria. All background isotopes are measured with the single electron (le) and electron plus 
gamma (ey) channels. A global analysis of the external background is discussed in |12J. -^^^Tl and ^''^Bi were inde- 
pendently measured using eyy, ejjj, and ea channels. The so-called "external background model" has been tested 
and validated using the dedicated sectors of ultra-pure Cu and Te foils in NEMO-3. 

Limits have been placed on the internal background activities of ^^Zr by a high purity germanium (HPGe) detec- 
tor, but ultimately the internal background activities are measured with the NEMO-3 apparatus. Internal background 
activities are measured in the le and ey channels. The le selection criteria are the following: one negatively charged 
particle track with length greater than 50 cm originating from the ''^Zr source foil and terminating at a scintillator, and 
an energy deposit > 500 ke V in the scintillator associated with the track. The ey selection criteria are the following: 
one negatively charged particle track with length greater than 50 cm originating from the '^^Zr source foil and termi- 
nating at a scintillator, an energy deposit > 200 keV in the scintillator associated with the track, an energy deposit 
> 200 ke V in a separate scintillator with no associated track, the cosine of the angle between the electron and gamma 
must be < 0.9 , and the time-of-flight information must be consistent with the electron and gamma originating from 
the same point in the source foil. In both the le and ey channels (for quality control of the reconstructed track) we 
require at least one triggered Geiger cell in first two layers closest to the source foil and less than 3 triggered Geiger 
cells that are not associated with the reconstructed track. 

The internal background activities are distinguished and measured due to contrasting Qp values and le and ey 
energy spectra of the isotopes. Equilibrium within a decay chain implies specific isotope activities to be correlated. 
^^^Ac, ^'^Bi, and ^os-pj ^.g p^.j jj^g 232'pjj chain and separated by short half-lives, therefore ^^**Ac and ^'^Bi activities 
are set equal and ^o^-pj jg ggj branching ratio of 36%. ^'^Bi and ^'"'Pb belong to the ^^^U chain and are set equal. 
234mp^ is also part of the ^-"^U decay chain but equilibrium with ^'"^Bi cannot be assumed due to the large half-life of 
the intermediate isotope ^^^Ra. Within this background model, contributions from the above isotopes to the le and 
ey channels have been fitted to experimental data over the entire energy region leaving the activities of the isotopes 
floating. Figure m shows the goodness of fit of the le channel and has a;t-^ = 93.9/45. Figure [3] shows the ey 
channel and has d^x^ - 43.7/28. The individual and summed energy distributions of electrons and gammas as weU 
as the angular distribution between them are plotted. The measurements of the internal ^^Zr contamination obtained 
in the le and ey channels compared with previously obtained HPGe limits in Table [T] provide a cross-check for the 
NEMO-3 measurements. The obtained numbers are in agreement with the ^'"^Bi and Tl activities (0.17 ±0.05 and 
0.08 ± 0.01 mBq respectively) reported in 1 12 1 where more restrictive energy intervals and diff'erent event topologies 
were used to identify signatures of the isotopes. 

Table 1 : Internal contamination of the '^Zr foil measured with NEMO-3 in the le and ey channels under the assumptions of the background model 
described in |5.1| The NEMO-3 measurements are compared to HPGe limits at 95% confidence level. 



The adjacent jSjS source isotopes ( Nd and Ca) and their associated internal backgrounds contribute events 
that pass the ^^Zr selection criteria due to the positional resolution of the Geiger cell tracker and accuracy of the 
reconstructed event vertex. These events have been studied and contribute ~ 1% in the le channel and ~ 7% in the ey 
channel and are included in the background description for ^^Zr 



Isotope 



NEMO-3 (mBq) HPGe (mBq) 




0.25 + 0.02 <0.75 

0.25 + 0.02 <0.75 

0.091 + 0.007 <0.23 

0.19 + 0.02 <0.45 

0.19 + 0.02 <0.45 

19.7 + 0.1 <19 

0.49 + 0.01 <6.6 
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Figure 2: Energy spectra of the Zr backgrounds in the le channel. Individual internal backgrounds are plotted (a) and the total background (b) is 
divided into 2 sub-groups of summed internal (int) and external (ext) components. 



6. Results 

6.1. Ivpp Decay 

The selection criteria for ee events are the following: two negatively charged particle tracks with lengths greater 
than 30 cm, both tracks originating from the '^Zr foil and terminating at independent scintillators, energy deposits 
> 200 keV in the scintillators associated with the tracks, each track has at least one triggered Geiger cell in first 
two layers closest to the source foil, and the time-of-flight information must be consistent with the two electrons 
originating from the same point on the source foil. 

The distributions of the energy sum of the two electrons, energies of the individual electrons, and the angle between 
two electrons are shown in Figure |4] 898 data events have been selected after 1221 days of data taking with a total 
expected background of 437.6 + 7.2 events. A maximized binned log-likelihood fit to the energy sum spectrum is 
performed to estimate the Ivfip signal contribution. The likelihood fit predicts 429.2 + 26.2 signal events (signal- 
to-background of 0.98) with a 7.5% efficiency. The breakdown of individual background contributions is shown in 
Table 121 

6.2. The Systematic Error 

The systematic error on the 2vy6y6 measurement has been investigated. The main contribution is from the error 
on the tracking detector resolution and track reconstruction efficiency |7 |. There is a 2% uncertainty in the mass of 
^*Zr fT\. The precision of the energy calibration of the calorimeter is 1% and the effect was determined by coherently 
changing the gain of the PMTs ± 1 % and observing the change in half-life. 

The systematic uncertainty of the external background model is considered. ^'"^Bi and ^os-pj -^^ j-j^g tracking 
chamber show a discrepancy between the channels they are measured in. ^'"^Bi is measured in the ey and ea channels 
and the obtained values differ by ~ 10% lfT2l . 208rpj measured in the eyy and eyyy channels and the obtained values 
differ by ~ 10% |12|. A conservative estimation on the total external background uncertainty is therefore 10% and 
is evaluated by fluctuating the external background + 10%. The attributed uncertainty on the measured half-life is 
±0.3%. 

The systematic error on internal background model is estimated by the difference in measured activities in the le 
and ey channels. The difference never exceeds 5% for the internal backgrounds, therefore the uncertainty on the 2vy6y6 
half-life is estimated by fluctuating the internal backgrounds +5% and recording the corresponding change in Ivfifi 
half-hfe. 

5 
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Figure 3: The ey channel displaying (a,b) the summed energy Eg + Ey, (c) the angular distribution between the electron and gamma cos(9), and 
(d) the energy of the electron Eg. As in Figure[2] the background contributions are divided into 2 sub-groups of summed internal (int) and external 
(ext) components. 



The world's best 2vySy6 half-life measurements for '^"Nd lfT6l fTTIl and '^^Ca ifTT l have been recently obtained. 
These isotopes neighbor the ^^Zr source and are included as backgrounds. The uncertainty in their measured half- 
lives is applied and the change in the ^''Zr half-Ufe is noted. The 2vyS/3 half-life of '^"Nd is known to 10% (including 
statistical and systematic errors) and contributes a ± 0.7% error on the obtained ^''Zr half-life. The 2vy6y6 half-life of 
*^Ca is known to 18% (including statistical and systematic errors) and contributes a + 0% error on the obtained ^^Zr 
half-life. 

"'"K is the dominant background in the ee channel and a systematic effect is observed by changing the energy 
window of the likelihood fit to exclude energy sums below 1.1 MeV. The strict energy window suppresses "^"K events 
and reduces the half-life dependence on the activity of ""^"K. The obtained systematic uncertainties are listed in Tab.js] 
and give a total systematic error of +6.7% and -6.2%. The final result for the 2v^/3 half-Hfe of ''^Zr including 
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Table 2: The number of events expected for the Zr internal and external backgrounds in the ee channel for 1221 days of runtime. 



Background 


Expected Nhkg 


Eff. (%) 




11.1±0.9 


0.042 


212Bj 


9.6 ±0.7 


0.036 


208 rpj 


9.3 + 0.7 


0.098 


214Bj 


22.8 + 2.5 


0.12 


214pb 


3.3 + 0.4 


0.017 


40k 


280.0 + 2.4 


0.014 




38.3 ±0.7 


0.074 


^^Ca internals 


0.0 + 0.0 




'^°Nd internals 


37.6 + 3.2 




External 


25.6 + 5.2 




Total 


437.6 ± 7.2 





Table 3: Summary of systematic errors pertaining to the 2y/Jy3 measurement of Zr. 



Description 


Syst. Error (%) 




the tracker and reconstruction 


±5.0 


|7| 


+ 1% energy calibration precision 


+ 2.9,-2.2 




the mass of ^^Zr 


±2.0 


m 


±10% external background precision 


±0.3 


mi 


±10% '5°Nd precision 


±0.7 




±5% internal background precision 


±1.9 




the likelihood fit energy window 


+ 1.6,-0.2 




Total Systematic Error 


+ 6.7%, -6.2% 





statistical and systematic errors is 

= [2-35 ± O.lMstat) + 0.l6(syst)] x lO'^r. (2) 

For comparison, (j2| is consistent and ~4 times more precise than the previous direct measurement (2.l^Q^(stat) + 
O.Ksyst)) X 10''' yr L18J . 

6.3. OvyS/3 Decay ( neutrino mass mechanism) 

Limits on Ov processes are evaluated using a binned log-likelihood ratio (LLR) test statistic lfT9l . Signal-plus- 
background and background-only hypotheses are defined and populated with simulated outcomes via Poisson statis- 
tics. The confidence level on the signal (CL,) is defined as CL, = CLs+h/CLt, , where CLs+i, and CL/, are the con- 
fidence levels of the signal-plus-background and background-only hypotheses respectively. The limit is found by 
scaling the signal until CLs reaches the 90% level. The detection efficiency for observing the Qv/3f5 decay is 20% and 
the limit at the 90% CL (as seen in Fig. |5(a)[ ) is 

rfjj >9.2xl02'yr. (3) 
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Figure 4: The (a) energy sum of both electrons E[ + E2, (b) MC signal fit to background subtracted data, (c) angular distribution cos((9), and (d) 
individual electron energy for 1221 days of runtime in the ee channel. The data are described by the sum of the expected backgrounds from MC 
and the 2v/J/J signal from the maximized log-likelihood fit. 



6.4. Other OvySjS Transitions 

More exotic mechanisms include Ov excited states, right-handed currents (RHC) and the emission of one or two 
Majoron particles (x'-'). The term "Majoron" is generally used to describe massless or light bosons which weakly 
couple to neutrinos. The various Majoron models are referred to by their spectral index n which determines the 
dependence of the phase-space on the energy released in the decay as G"'' oc (g/3/3 - T)", where T is the sum of the 
kinetic energies of the electrons emitted in the decay. The standard model 2v/3/3 decay has n - 5. For single Majoron 
emission (x°/3/3) n = 1 (201 |2l|, for double Majoron emission CrVySyS) n = 3 or n = 7, 1221 123] |23 and for the 
"bulk" Majoron in the context of the "brane-bulk" scenario n - 2 Il24ll . The limits have been obtained using the 
binned log-likelihood ratio test statistic method described in 6.3 The results are presented in Tables]?] and |5] Naively 
one would expect the mass mechanism and RHC limit to be identical because they have the same summed energy 
spectrum, but due to the asymmetrical energies and small angle (cos 9 ~ \) between the emitted electrons in the RHC 
process, detector efficiencies decrease the RHC acceptance. 




Ei + Ej (MeV) Ei + Ej (MeV) 



(a) OvjS/S mass mechanism (b) « = 1 Majoron emission 

Figure 5: Energy sum distributions of two electrons Ei + E2 for 1221 days of runtime. The data are described by the sum of the expected 
backgrounds from MC. the measured 2v/3/3 signal, and the limit on (a) the OvjSjS mass mechanism and (b) the n = I Majoron emission at 90% CL 
using the log-likelihood ratio test statistic method described in |6.3| The Ov histogram is not included in the total MC (dark-grey) 

Table 4: Summary of half-life limits T1/2 (yr.) evaluated at the 90% CL for Ovflfi mechanisms where O^^siimpp)) is the standard Ov/?/3 decay to the 
ground state, d^dmpp)) is the first excited state, Ojj(RHC;/l) is the right-handed cuiTent decay to ground state and 2j^(RHC;/l) is the first excited 
state. The bottom line is the previous best limit presented for comparison. 



OvyS/3 






0;,(RHC;^) 


2|(RHC;/1) 


9.2 X lO^i 


2.2 X lO^" 


5.1 X lO^i 


9.1 X lO^o 


i.ox lo^'insi 


6.8 X IQi'* 


Ea - 


3.9 X 10^°JT8ll 



Table 5: Summary of half-life limits 7'i/2 (yr.) evaluated at the 90% CL for Majoron emission decay mechanisms. The spectral index («) for the 
Majoron modes refers to the dependence of G"'' oc (Qpp — T)" where T is the electrons' kinetic energy sum. The bottom line is the previous best 
limit presented for comparison. 



Majoron 


n = 1 


n - 2 n — 3 


n = 7 


1.9 X 10^1 


9.9 X 10^0 5.8 X lO^o 


1.1 X 10^° 


3.5 X lO^OES 


1 - 6.3 X lO^^ESl 


5.1 X IQi'^lilSI 



7. Discussion 

7.7. IvNME 

The largest uncertainty in the effective Majorana mass determination is due to the uncertainty of the Ov/3/3 NME 
(M^^). It is still difficult to calculate the NMEs for '"'Zr and currently there are no large-scale shell model calcu- 
lations (see review ll27l ). The current models for the M"'' computation of ^^Zr are the quasi-particle random-phase 
approximation (QRPA) and the renormalized (RQRPA) |4 5 6 1, but unfortunately cannot precisely predict M"'' due to 
strong dependence on the unknown parameter gpj, (particle-particle coupling). In fact, extracted experimental values 
of M^" are needed to fix gpp which is used for the Af"" computations. Two values for the parameter gy^ are generally 
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agreed upon and the NME is computed using both gA — 1.0 and gA — 1.25. Recently a new approach (Projected 
Hartree-Fock-Bogoliubov - PHFB model) was developed ll28l l29l l30l which can predict the M^^ and M"'' values. 

Using the measured value of the ^^Zr 2v/3/3 half-life ^ we extract the experimental value of the corresponding 
NME according to the formula 

[T^J^]-^ = G^^IM^f , (4) 

where G^^ = 1.8 X 10"^^ yr"' is the known phase-space factor ll27l using gA = 1.25. The extracted NME is scaled by 
the electron rest mass and is 

M^" = 0.049 + 0.002 . (5) 

One can compare this (jsjl value with the calculated value, M^" = 0.058 128 J . The obtained precise value for M^'^ will 
be used to fix gpp parameter and improve the M^'' calculations for ^*Zr 

7.2. Majorana Neutrino Mass 

The limit on the Qv/3/3 half-life ([3]) is used to calculate an upper bound on the effective Majorana neutrino mass, 
{nipp} from Using the NME computed in QRPA f^l and (R)QRPA |5 6| the following limits on the eff'ective 
Majorana neutrino mass are obtained. The range is due to the M"'' uncertainty: 

(mpp) < 7.2 - 8.2 eV (NME from gl) , (6) 
{t7ip/s) < 11.0 - 19.5 eV (NME from BISl) . (7) 

Using the NME from |29 | (PHFB approach) we obtain (m^j^) < 15.1 eV. 

7.3. Majoron-Neutrino Coupling 

A limit on the Majoron to neutrino coupling {gee) is extracted for the most commonly discussed mechanism with 



n = 1 (as seen in Fig. 5(b) i from 

[r°f]-'=G'"'^°|M»>'^?<g,,)2. (8) 

The Majoron matrix element M"'*" is identical to that of the mass mechanism (M"") and using NME from QRPA ID 
and (R)QRPA ||5] |6l the following limits on the Majoron coupling are obtained. The range is due to the M^^-^" 
uncertainty: 

<g,,)<(1.5-2.1)xl0"'* (NME from 11), (9) 
<g,,) < (3.3 - 5.7) X lO""* (NME from [Dig). (10) 

Using the NME from ||29ll (PHFB approach) we obtain {gee) < 4.4 x 10 For mechanisms n = 2, 3, and 7 our limits 
are comparable with the best present limits obtained with other nuclei IJTI . 

7.4. Gf Time Variation Hypothesis 

It has been suggested in l32l [33l that observed differences in half-lives of fifi isotopes obtained in geochemical 
experiments with samples of different age could be related to time dependence of the Fermi constant Gp. Due to the 
stronger dependence on the Fermi constant (G^ rather than G^), decay offers a better sensitivity than single (3 decay 
studies. The ^^Zr-'''Mo transition is of particular interest since the daughter element is not a gas thus eliminating the 
main systematic error of the geochemical measurements. A comparison between the half-lives obtained with ancient 
zircon (ZrSi04) minerals characterizing the decay rate in the past with present day fifi decay rates obtained in a direct 
experiment like the one presented here allows the hypothesis to be probed with a high precision. 

A previous geochemical measurement carried out in 1992 with a 1.7 x 10^ yr old zircon yielded a 2vyS/3 half-life 
of (3.9 + 0.9) X lO'^ yr l34l . An independent measurement was performed in 2001 with a number of zircons aged 
~ 1.8 X 10^ yr and a half-life of (0.94 ± 0.32) x 10'^ yr was measured |35|. The measurement presented in this paper (j2]i 
lies between the two geochemical measurements. More accurate studies of minerals of different age are needed in 
order to probe the Gf time variation hypothesis with high precision. 
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8. Summary 



The most precise measurement of the 2vj6y6 decay half-Hfe of Zr to date has been presented including the charac- 
teristics of the final state electrons (energy sum, individual electron energy, and angular distribution). The presented 
limit on the Qvfifi half-life Q is an order of magnitude higher compared to the previous result 1 18|. The presented 
limits on other Ov processes are the best to date. The high Qpp value, hence large phase-space of ^^Zr, makes it an 
excellent choice for the study of OvySyS decay if the enrichment of this isotope in large quantities proves to be feasible. 
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